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Abstract 

We study the equilibrium domain structure and magnetic flux around a 
ferromagnetic (FM) film with perpendicular magnetization Mq on a super- 
conducting (SC) substrate. At 4ttMq < H c \ the SC is in the Meissner 
state and the equilibrium domain width in the film, I, scales as (I/AttXl) = 
(In /4:7tXl) 2 ^ 3 with the domain width on a normal (non-superconducting) 
substrate, In/4tt\l 3> 1. Here Xl is the London penetration length. For 
4ttMq > H c i and In in excess of about 35 Xl, the domains are connected 
by SC vortices. We argue that pinning of vortices by magnetic domains in 
FM/SC multilayers can provide high critical currents. 



The interaction between magnetism and superconductivity has been intensively studied 
in the past, see, e.g., the revie J. The discovery of high-temperature superconductors and 
advances in manufacturing of nanoscale multilayered systems have added new dimensions to 
these studies. In this paper we investigate equilibrium magnetic and superconducting phases 
in a system consisting of a ferromagnetic (FM) film with perpendicular magnetic anisotropyi 
on the surface of a superconductor (SC). The interest to such systems is two- fold. Firstly, 
the needs of magneto-optic technology have produced a large variety of magnetic films with 
perpendicular magnetic anisotropy. Some are synthesized on metallic substrates (e.g., Nb) 
and are well-characterized at room temperature, including their domain structurei'SII. We 
will show that, as the temperature of such a system is lowered, its magnetic state must 
be affected, in a non-trivial manner, by the superconducting transition. Secondly, in a 
multilayered SC/FM system, the domain structure in the FM layers should produce the 
pinning of superconducting vortices which may be significantly stronger than the pinning 
by magnetic dotsi0ii. 

The system under consideration is shown in Fig.l. We are assuming no exchange of 
electrons between the FM film and the superconductor. This will be true either when 
the ferromagnet is an insulator or when it is separated from the superconductor by a thin 
insulating buffer layer. Then the FM film and the superconductor are coupled only by the 
magnetic field. (The systems with the exchange of electrons between the FM and SC layers 
have been discussed in Ref. |I0| .) In this case the superconductivity makes a profound effect 
on the domain structure in the FM layer, which must be easy to detect in experiment. The 
physics behind this effect is explained below. Consider a FM film of thickness (1m, with 
perpendicular magnetic anisotropy. In the absence of the superconductor adjacent to the 
film, its domain structure is determined by the balance of the energy of the magnetic field 
surrounding the film and the energy of domain walls. The positive energy of the magnetic 
field favors small domains, so that the field does not spread too far from the film. On the 
contrary, the positive energy of domain walls favors less walls, that is, large domains. The 
minimization of the total magnetic energy gives a well-known result0 for the equilibrium 



domain width, Ioc^/Mm, with 5 being the domain wall thickness. Domains typically observed 
in magneto-optic films have thickness of a few micron. In the presence of a superconductor 
adjacent to the FM film, the balance of the magnetic energy changes drastically. This is 
because the magnetic field must be either expelled from the superconductor due to the 
Meissner effect or it should penetrate into the superconductor in the form of vortices. In the 
first case, the superconductor favors FM domains of width below the London penetration 
depth, Al- If the room temperature domains are significantly greater then \l, the effect of 
the SC phase transition on the domain structure will be dramatic. As we shall see, the new 
equilibrium will be achieved at 1(x(5(Im^l) 1 ^ '■ Consequently, on lowering the temperature 
below the SC critical temperature, the domains in the FM film can shrink by an appreciable 
factor. 

We are assuming the stripe domain structure in the FM film. The width of the FM 
domain, /, is presumed large compared with the domain wall thickness 5. The latter is the 
smallest length in our consideration. Two other characteristic length are the thickness of 
the FM film, d,M, and the London penetration depth, Al, of the SC. In the case of / < A^ 
the magnetic flux penetrates into the SC as it would penetrate into a normal non-magnetic 
metal, making superconductivity irrelevant. The case of interest is, therefore, / > \ L . We 
shall begin with the study of the Meissner state, that is, the state where equilibrium vortices 
are absent. 

The free energy functional for the magnetic field, B = [B x (x, z), 0, B z (x, z)}, is 



Here J-'s is the free energy due to the magnetic field in the superconductor, Tm is the free 
energy of the magnetic film and the empty space above the film, M(x) is the magnetization 



^(B,M) = ^ 5 (B,M)+^ M (B,M) , 



(1) 



where 



^5(B) = ^/ ^[A|(VxB) 2 + B 2 ] 
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inside the magnetic film, and Tb is the energy of domain walls. At / >> 5, a good approxi- 
mation for M(x) (see Fig.l) is the step-like function along the X-axis, M(x) = ±M inside 
the domains. Its Fourier expansion is 

where M is the magnetization and Q = 2n/l. For this domain structure !Fb(1) = ctcIm/I- 
Here a = V2/3Mg6/ 7r is the energy of the unit area of the domain wall and (3Mq is the 
energy density of the perpendicular magnetic anisotropy. 

The equilibrium distribution of the magnetic field should be obtained by the minimization 
of J^(B, M) at a given configurations of magnetic domains M(x). Introducing H = B— 47rM, 
one obtains in terms of H: 



^a/(H,M) 
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+ F D , Js(H) = F S (B = H) . (4) 



The field H is induced by alternating magnetic charges, V • M, on the two surfaces of the 
frnetic filirS With account of the Maxwell equation, V-B = 0, it satisfies 



ma 



V • H = -4ttV • M = -An[6(z) - 5{z + d M )\M{x) , V x H = (5) 

outside the superconductor, that is, inside the magnetic film, in the buffer layer, and in 
the empty space above the film. Here 8(z) is the delta-function, z = and z = —du are 
coordinates of the film surfaces. Inside the superconductor H satisfies the London equation, 

V 2 H-AzT 2 H = , (6) 

and the boundary condition that H is continuous across the interface between the buffer 
layer and the superconductor. Equation (6) is valid if the magnetic field changes on the 
scale greater than the correlation length £. In our case the smallest relevant scale of spatial 
variations of the field is the width of FM domains I. We shall assume that I » £, which is 
relevant to most situations of practical interest. 
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Solving the above equations we obtain that, due to the domain structure, H(x, z) decays 
exponentially away from the surfaces. Taking into account that du S> I and neglecting 
exponentially small terms of order oc exp[— du (47r 2 Z -2 + A^ 2 ) 1 / 2 ], we get 

H(ar, z) = ^2U q exp(-q z \z\ + iqx) (7) 
q 

for the magnetic field inside the superconductor, the film, and in the empty space. Here 
z is the distance along the z-axis from the nearest film surface. This gives the Fourier 
components 

00 AM 

H z , g = j:^ T S[ q -Q(2k + l)} 

00 AM 

H z ,- q = ~ £ ^7% + Q( 2k + !)1 > ( 8 ) 

and H x>q = —iq z H z>q /q with q 2 z = q 2 + Xj 2 inside the superconductor and q z = q elsewhere. 
Substituting this equilibrium magnetic field H at a given I into the free energy functional, 
Eq. (f|), we obtain the following expressions for !Fs(l) an d ^Fm(1) per unit area: 

T(n 4M 2 - [1 + (2A;+1) 2 Q 2 A L 2 ] 1/2 

Fm(1) = 3Fs(l, AZ 1 = 0) + F D (l). (10) 

Above T c the free energy of the system as a function of I is given by 

T N (l) = AF s {l,\ L l =0)+F D . (11) 

The minimization of (11) gives the well known result for the equilibrium width of the domains 
when the superconductor is in the normal state0, 

nl/2 



N 



V2n 



(P5d M )V 2 . (12) 



_7C(3)_ 

For the superconducting state of the substrate it is convenient to introduce I = 1/AttXl 
and In = In/^^l- Then 



The minimization of T with respect to I produces the dependence of I on In shown in Fig. 2. 
At In <C 1 the field penetrates into the SC same way as it penetrates into the normal 
metal and Z ~ Zat (see insert to Fig. 2). In the opposite case of In 1, a rather accurate 
approximation is 

l^fjl 3 . (14) 

We, therefore, conclude that the SC phase transition in the substrate can result in a signif- 
icant shrinkage of the equilibrium domain width in the FM film if the substrate is in the 
Meissner state. 

The Meissner state studied above should always be the case when AttM < H cl . At 
47rM > H c i the equilibrium energy of the Meissner state, JF M , should be compared with the 
energy of the vortex state. Vortex lines connecting domains with opposite magnetization can 
form only if l N 3> 1. The energy of the SC in the vortex state may be easily estimated in the 
limit of M 3> 4:irH cl . In that case the average distance between vortices is small compared 
to their magnetic radius \l and the average field in the SC is close to that in a normal metal. 
The corresponding total free energy of the system is then close to jFjy. Some small corrections 
to that energy arise from the vortex line tension and from the repulsion of vortices, which 
is of order &/(<I>oMo/87rA|) ln(if c2 /47rM ) accounting for the curvature of vortex lines 

and for their additional energy near the FM surface). These corrections are small in the 
limit of large magnetization. Consequently, for the vortex state of the SC substrate, the 
equilibrium domain width in the FM film should be close to that on the normal substrate. 
For the ratio of the equilibrium free energies one obtains T )Tn = [G/7((3)]l N ' 3 m0.7l3l N ' 3 , 
where T is computed at I — f^ 3 and J-'n is computed at I — In- Thus, at 4ttM > H& and 
Tjv>2.8 (that is for l N in excess of about 35A^) the vortex state should be energetically more 
favorable than the Meissner state. 

It should be emphasized that all conditions, obtained in this paper, depend on tem- 
perature through the temperature dependence of A^ and In- Nevertheless, the factor that 
determines the maximum shrinkage of equilibrium FM domains in the Meissner phase is 
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a universal number, max(ZAr/Z)=max(r]v 3 )=7C(3)/6~1.4. This should be easy to detect in 
experiment. 

The effects described above fall within common experimental range of parameters. They 
will be noticeable if the room-temperature domains in the FM film are wider than Zjv~0.5 
micron for the Nb substrate (Al~40 nm) or wider than Ztv~1.6 micron for a high-temperature 
SC (A^~130 nm). Because equilibrium l N depends on the thickness of the FM film, o? M , the 
above condition on translates, through Eq. (12), into the lower bound on du- For, e.g., 
a TbFe film! (/?~10 2 and <5~15 nm) the equilibrium domain width should decrease below 
T c in films of thickness greater than 0.3 micron on a Nb substrate or in films thicker than 3 
micron on a high-T c substrate. For TbFe film on a Nb substrate the Meissner state should 
occur for l N <1.4 micron (<i M <2.4 micron) and the vortex state should occur at 1^ >1.4 
micron (d,M >2.4 micron). In the case of a high-T c substrate the Meissner state should occur 
for 1^ <4.5 micron (d,M <25 micron), while at In >4.5 micron (g?m >25 micron) the vortex 
state should occur. 

In a real FM film the stripe domains are curved due to thermal fluctuations^ and due to 
the pinning of domain walls. This, however, should not affect our conclusions as long as the 
corresponding radius of curvature of domains is large compared with other characteristic 
lengths. Since we are interested in the equilibrium magnetic structure due to the FM- 
SC interactions, it is important to acknowledge that strong pinning of domain walls by 
the imperfections may prevent the system from reaching that equilibrium. Possible ways to 
study the equilibrium magnetic structures include choosing systems with low coercivity (that 
is, weak pinning of domain walls), or low Curie temperature (below the critical temperature 
of the SC), or rotating the system in a slowly decaying magnetic field. It should be also 
possible to extract changes in the magnetic equilibrium from the study of the magnetic 
hysteresis in the FM film above and below T c . A large variety of magnetic materials should 
allow experiments in all interesting ranges of temperature and coercivity. 

Finally, we would like to comment on the magnetic pinning of vortices in SC/FM mul- 
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tilayers. In the past the enhancement of pinning in conventional superconductors was done 
through manufacturing samples with microscopic defects. Several new approaches have been 
developed in recent years. They include manufacturing films with micrometer-size holes0 
and depositing magnetic particles^! or magnetic dotsS onto SC films. For high temperature 
superconductors a remarkable effect has been achieved in samples with columnar defects 
produced by heavy ion irradiation^. All the above methods are based upon introducing 
defects that suppress superconductivity. The pinning then arises from the tendency of the 
vortex normal core to match with the region where the superconductivity is suppressed. The 
maximal energy per unit length of the vortex, available for such pinning, is the condensation 
energy of Cooper pairs in the volume of the vortex core, (if ( ?/87r)7r^ 2 ^($ /8 7r ^L) 2 - It is easy 
to see that the pinning of the magnetic flux of vortices in SC/FM multilayers can be signif- 
icantly stronger than the pinning of vortex cores. As the problem of statistical mechanics 
it will be reported elsewhere@. Here we shall just estimate the amplitude of the one-vortex 
pinning potential due to magnetic domains. This estimate can be obtained by considering 
the energy of an extra vortex, created by an external magnetic field in the presence of the 
domain structure in FM layers. The upper bound on the magnetic pinning energy is <&qMq 
per unit length of the vortex. For Mo=500 emu/cm 3 it is one hundred times the energy of 
the pinning by columnar defects. The magnetic pinning of SC vortices will be effective if 
the field in the FM layers does not exceed the coercive field. Above that field the pinning of 
domain walls disappears and FM domains, together with SC vortices, become mobile. Thus, 
strong pinning of vortices favors large coercivity of the FM film, the condition opposite to 
the one required to observe the shrinkage of FM domains in the Meissner phase. Some 
evidence that magnetic dots produce more pronounced pinning than non-magnetic dots has 
been recently demonstrated in Nb films@. Another study0 indicated that the irreversibility 
line of an YBCO film is pushed up when a barium ferrite film is deposited on the surface 
of the SC film. Thus, the study of SC/FM multilayers, besides being a fascinating "mul- 
tidimensional" problem, can also be promising in developing SC systems with high critical 
currents. 



This work has been supported by the U.S. Department of Energy through Grant No. 
DE-FG02-93ER45487. 
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FIGURE CAPTIONS 
Fig. 1. 

FM film with stripe domains on a SC substrate. 
Fig. 2. 

Dependence of I (normalized equilibrium domain width on the SC substrate) on l N (nor- 
malized domain width on the normal substrate). 



12 




z 



